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It is commonly accepted that there are no phase 
transitions in one-dimensional (ID) systems at 
a finite temperature, because long-range correla- 
tions are destroyed by thermal fluctuations. Here 
we demonstrate that the ID gas of short-range in- 
teracting bosons in the presence of disorder can 
undergo a finite temperature phase transition be- 
tween two distinct states: fluid and insulator. 
None of these states has long-range spatial corre- 
lations, but this is a true albeit non-conventional 
phase transition because transport properties are 
singular at the transition point. In the fluid phase 
the mass transport is possible, whereas in the in- 
sulator phase it is completely blocked even at fi- 
nite temperatures. We thus reveal how the inter- 
action between disordered bosons influences their 
Anderson localization. This key question, first 
raised for electrons in solids, is now crucial for 
the studies of atomic bosons where recent exper- 
iments have demonstrated Anderson localization 
in expanding very dilute quasi-lD clouds. 

The absence of finite temperature phase transitions in 
one-dimensional (ID) systems is considered as almost a 
dogma. Its justification is based on another dogma which 
states that any phase transition is related to the appear- 
ance/disappearance of a long-range order or at least long- 
range spatial correlations. Thermal fluctuations in ID 
systems destroy any long-range order, lead to exponen- 
tial decay of all spatial correlation functions and thus 
make phase transitions impossible^. Non-interacting 
quantum particles in a one-dimensional random poten- 
tial show a similar behavior in the sense that all single- 
particle eigenfunctions are localized, i.e. decay exponen- 
tially in spaced The same statement holds for two di- 
mensional systems without spin-orbit interactions^. By 
contrast the single-particle states in three dimensions are 
either localized or extended as a result of the Anderson 
transition^. In this paper we demonstrate that the ID 
gas of weakly short-range interacting bosons in the pres- 
ence of disorder exhibits a finite temperature phase tran- 
sition between two distinct states, fluid and insulator, 
and the transition temperature depends on the disorder. 
None of these states is characterized by a long-range or- 
der or long-range spatial correlations. Moreover, ther- 
modynamic functions, such as specific heat, do not have 
singularities at the transition point. From this point of 
view, the dogma is not violated. Nevertheless, this is 
a true albeit non-conventional phase transition, because 
transport and energy dissipation properties of the fluid 



and insulator phases are dramatically different and are 
singular at the transition. The difference between the 
fluid and insulator phases can be qualitatively under- 
stood by comparing two many-body ID systems: inter- 
acting particles without disorder and the ID Anderson 
insulator (disorder without interactions). In the fluid 
without disorder, the dissipation of energy of an arbi- 
trarily slow external field and mass transport are pos- 
sible. At the same time, even at finite temperatures 
there is no mass transport in Anderson insulators and 
the energy dissipation vanishes for the frequency of the 
external field tending to zero. Here we show that inter- 
acting ID bosons in disorder demonstrate one of the two 
types of behavior and describe the phase diagram in the 
temperature-disorder plane. We thus provide an answer 
to the subtle question of how the interaction between 
disordered particles may suppress Anderson localization 
and permit them to acquire the fluid behavior. This 
was the key problem for charge transport in electronic 
systems, and it is now emerging in a new light in the 
studies of disordered ultracold bosons. These studies are 
driven by fundamental interest and by potential applica- 
tions of atom waveguides on a chip^. Recent remarkable 
experiments^^ have demonstrated Anderson localization 
in expanding extremely dilute quasi- ID Bose gases, and 
the investigations of effects of the interparticle interac- 
tion in relatively dense clouds are underway. 



I. SINGLE PARTICLE LOCALIZATION IN ID 

Let us first discuss the density of states (DoS) for a 
single particle with mass m in ID and introduce relevant 
energy and distance scales. In the absence of disorder the 
particle eigenstates are plane waves with energies e > 0. 

The DoS is Va(e) = \J m/2ir 2 h 2 e and it diverges at e — > 0. 
What happens in a static random potential? At large and 
positive e the DoS is only slightly affected by disorder: 
v(e) w vo{e) ~ 1/y/e. At the same time, the disorder 
cuts off the DoS divergence at e = and transforms it 
into a peak with a finite height and width E„ as shown 
in Fig. [I] Also, the disorder creates states with negative 
energies, thus making the DoS finite at e < 0. For large 
negative e the DoS is exponentially smal l, and this part 
of the curve v(e) is known as Lifshitz tai l 9 ' 10 * n 4 

Both v„ and E* are determined by the statistics of 
the random potential U{x). For simplicity we assume 
a short-range Gaussian potential with the amplitude 
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Uq and correlation length a such that condition Uq <C 
ft 2 /ma 2 holdJ ^. The n the only relevant energy and 
length scales ar d 9 * 10 * 11 -!: 

E, ~ (Ctfo'm/ft 2 ) 1 / 3 ; (1) 
C* = h/y/m~K = (tf/Ugam 2 ) 1 / 3 . (2) 

They determine the width of the DoS peak and the 
maximum DoS value i/„ ~ 1/E*(*. In order to obtain 
these scales consider a weakly bound state of a parti- 
cle in the potential U(x), with an extension of the wave 
function, £ ^> c The particle energy can be written as 
E ~ (h 2 /2m( 2 — Uoy/a/Q, where for the Gaussian dis- 
order the potential energy term is obtained multiplying 
the contribution of each potential well, Uqu/C,, by \J (/a, 



which is the square root of the number of wells on the 
length scale £. The energy E reaches a minimum value 
E* at C ~ C* 5 w ith E„ and C* given by Eqs. (Q and 

We thus see that the single-particle spectrum can be 
divided into three parts: high-energy states with e E* 

and v(e) ~ \J m/h 2 e, low-energy states located in the 
region of the DoS peak, and Lifshitz tail at negative en- 
ergies. As we already noted, all single-particle eigenfunc- 
tions in ID are localized with an energy-dependent lo- 
calization length ((e). For high-energy states, e 3> E*, 
we have £(e) ~ e(*/E*, whereas in the Lifshitz tail 
C(e) ~ (*\/Ej\ef^^. For the low-energy states in the 
DoS peak, |e| < E„, the localization length is £(e) ~ £*■ 




FIG. 1: Properties of single particle localization. In (a) the brown curve shows the disorder potential U(x), with blue solid lines 
and dashed red lines indicating the location of tail and high-energy states, respectively. The shape of the wavefunctions of these 
states is shown by blue (tail) and red (high-energy) curves tp(x). In (b) the density of states v and localization length £ versus 
energy e are shown in blue for the DoS peak and tail states, and in red for high-energy states. For studying the many-body 
localization transition, v(e) and ((e) in the DoS peak (low-energy states) can be approximated by the green dotted lines. 
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II. MANY-BODY 
LOCALIZATION-DELOCALIZATION 
TRANSITION 



Repulsive short-range interaction between bosons in 
ID gives rise to two other energy scales in addition to E* : 
the temperature of quantum degeneracy = fi 2 n 2 /m 
(we use the units with Boltzmann constant ks — 1) , and 
the interaction energy per particle, ng, where g is the 
coupling constant for the interaction, and n is the mean 
boson density. We focus on the weakly interacting regime 
where the dimensionless coupling strength is small: 



7 = ng/T d < 1 



(3) 



so that at the mean separation between the bosons their 
wavefunction is not influenced by the interactions. It is 
convenient to introduce the dimensionless temperature 



t = 



T 
ng 



1 T 



(4) 



Another dimensionless parameter characterizes the 
strength of the disorder: 



E* I ng 



(5) 



so that large values of k correspond to strong disorder. 

At this point we should make two important state- 
ments. If the disorder is extremely strong (k — > oo), 
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the bosons occupy only states in the Lifshitz tail. As 
the DoS in the tail is exponentially small, the bosons 
are distributed among the "lakes" located exponentially 
far from each other. The bosons can not hop between 
the lakes, and the system is in an insulating state. The 
reduction of n at T — eventually (for k, — k c ~ 1) trans- 
forms the insulator into an algebraic superfiuid (spatial 
phase correlations do decay, but only algebraically) . This 
Kosterlitz-Thouless type transition was first analyzed in 
RefP and more recently discussed in relation to disor- 
dered Josephson chains^ and to cold atomic gase J 15 * 16 - 1 . 
However, well before the transition from insulating to 
superfiuid state most of the particles find themselves in 
low-energy states where the DoS is much higher than in 



the Lifshitz tail. Thus, we may neglect the tail in our 
discussion of the fluid-insulator transition and consider 
only low-energy and high-energy states (green and red in 

Fig. SO- 

Second, although the interaction between bosons 
renormalizes (screens) the disordered potential U(x), this 
does not change the picture of single-particle eigenstates 
in our discussion. The reason is that relevant particle 
energies are of the order of E* or larger. For k 1 
they exceed the interaction energy ng, and as we will see 
there is no need to consider k« 1. In this respect, the 
main effect of the interaction on the ground state of the 
system is not screening the random potential but rather 
controlling the occupation of single-particle states. 
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FIG. 2: Scattering processes leading to the many-body localization-delocalization transition for a classical Bose gas (T > Td) 
in (a), for a degenerate thermal Bose gas (Td^y < T < Td) in (b), and for the low-temperature regime (T < Td^/j) in (c). 



At a finite temperature T it is crucial to take into ac- 
count two-body processes that change occupation num- 
bers and can dramatically affect the properties of the 
system. For example, the system of interacting localized 
fermions (electrons) can have a finite DC conductivity 
even in the absence of coupling with any outside batrP^, 
whereas without interactions the conductivity is exactly 
zero at any T . In the presence of the interactions the con- 
ductivity remains zero unless the temperature exceeds a 
critical value T c . This transition can be thought oP^ 
as Anderson localization^ of many-body wavefunctions. 
The critical temperature T c depends on the interaction 
stre ngth. This is the many-body analog of the mobility 
edg d 18 * 19 -* that separates bands of localized and extended 
states in the single-particle Anderson transition. 

The Anderson transition is based on the fact that two 
quantum states belonging to different lattice sites hy- 
bridize provided that the hopping matrix element be- 



tween these states exceeds the difference in onsite en- 
ergies. As soon as the density of the hybridized pairs 
exceeds a critical value the eigenstates turn out to be 
extendecP. The many-body localization-delocalization 
transition can be qualitatively understood by extending 
this physical picture to states of more than one particle. 

Consider an occupied localized single-particle state \i) 
with energy e^. The interaction of a particle occupying 
this state with a particle in the state |fc) can transfer 
the |i)-state particle to the state \j), transferring simul- 
taneously the |fc)-state particle to another state \l) (see 
Fig. [2]). Let Uik.ji be the matrix element of this pro- 
cess. Due to the exponential decay of the localized wave 
functions one may assume that Uikji = unless all four 
sta tes a re located near each other. Moreover, it turns 
out! 20 * 21 ! that the matrix element rapidly decreases with 
an increase in the energy transfer |e$ — £fe|. We thus may 
confine ourselves to the case where the states \i) and 
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are nearest neighbors in the energy space. For simplic- 
ity we replace Uikji by a certain typical value Ut pro- 
vided that the states \i), \j), \k), \l) are localized nearby 
and pairwise are nearest neighbors in energy. 

In a random system the energies of the final and initial 
states can not be matched exactly. As long as the energy 
mismatch A^-ji — \ei+£k~ e j~ e i\ exceeds Uikji the effect 
of the interactions on the quantum state of 4 particles is 
negligible. One may say that single-particle excitations 
do not decajEH. Suppose that the interaction is weak 
and a typical mismatch A t exceeds the matrix element 
Ut. Does it mean that single-particle excitations have 
an infinite lifetime? The answer depends on the num- 
ber of channels, Ni, for the decay of a given excitation 
(more precisely, N\ is the number of possible processes 
\i), \k) — > \j), \l) that involve a given state \i)). Indeed, 
with probability of order unity, these processes should 
have a channel with the mismatch that is smaller by a fac- 
tor of Ni than the typical value A t . Therefore, Ut should 
be compared with A t /Ni. Note that N\ plays a role of 
the number of nearest neighbors in the single-particle lo- 
calization problem. However, in the many-body case N% 
is determined by the density of thermal excitations and 
is temperature dependent. Since characteristic single- 
particle energies and localization lengths are determined 
by the temperature, both A t and U t are also temperature 
dependent. As a result, there is a critical temperature T c 
following from the equation 

A t {T c ) = U t (T c )N 1 (T e ). (6) 

At T > T c many-body states are extended, i.e. they 
are linear combinations of one-, two-particle, etc. excita- 
tions, and the number of terms is infinite. This leads to 
the fluid behavior. For T < T c the many-body localiza- 
tion takes place and one expects the insulating behavior. 

Note that the arguments given above and Eq. (|6| are 
general and independent of quantum statistics of the par- 
ticles. At the same time, both Ut(T) and N\(T) do 
depend on the statistics. For disordered bosons these 
quantities are determined by the density of single parti- 
cle states v(e) and by the occupation number N(e). The 
latter is controlled by the chemical potential a which is 
related to the mean density n and temperature T by the 
normalization condition n = J deN(e)u(e), with N(e) de- 
termined by the Bose-Einstein distribution. More pre- 
cisely, as long as the interaction is weak (7 <C 1) the oc- 
cupation is N(e) = {exp[(e#F — n)/T] — where chf 
differs from the single-particle energy e by the Hartree- 
Fock corrections. We will see that these corrections be- 
come important only at sufficiently low T and are neg- 
ligible at high temperatures where in the vicinity of the 
transition relevant particle energies greatly exceed ng. 

III. PHASE DIAGRAM 



to represent the phase diagram in terms of the dimension- 
less temperature t = T /ng and find the dependence of t c 
on the strength of the disorder k. The relation t = t c (n) 
determines the boundary between the fluid and insulator 
phases in the (t, k) plane. Alternatively, one can speak 
of a temperature-dependent critical disorder n c (t). 

At high temperatures T > T4, or t > 7 -1 , the Bose 
gas is not degenerate and a characteristic energy of par- 
ticles is of the order of T ^> ng. As will be seen, near the 
transition temperature we have E* < T so that most 
of the particles are in the high-energy states and the 
occupation of all the states is small and described bu 
Boltzmann distribution. In this case the typical matrix 
element U t (T) in Eq. §6§ does not depend on the occu- 
pation and can be estimated as U t (T) ~ g/((T). The 
typical mismatch is the meen nearest neighbor energy 
spacing: A t ~ [^(T)C(T)]- 1 . The quantity Ni(T) is 
given by the number of particles localized within a dis- 
tance ~ C(T) from a given state \i), i.e. N^T) « ((T)n. 
Then, Eq. @ is reduced to ngv(T c )C,(T c ) ~ 1. Using the 
high-energy localization length £(T) = (*T/E* and the 
density of states v(T) = l/y/E^QT, with the help of 
Eqs. Q and (|5| we obtain the high temperature relation 
for the critical disorder : 

K c (t)~t 1 /8; t> 1 -\ (7) 

Close to the transition we have E*/T = n/t ~ i~ 2 / 3 <C 
7 2 / 3 <C 1 , which justifies our initial assumption that most 
of the particles are in the high-energy states. 

As the temperature is reduced below T d the Bose gas 
becomes degenerate. In the absence of disorder the chem- 
ical potential is u = —T 2 /Td as long as T > Td^, 
or t > 1/^/7. Characteristic energies of particles are 
~ \fi\ < T. However, as will be seen below, they still 
exceed both ng and E*. Therefore, most of the ID dis- 
ordered interacting bosons in the temperature interval 
7 -1 / 2 < t < 7 _1 occupy high-energy states and one 
can use the ideal Bose gas distribution N(e), i.e. neglect 
Hartree-Fock corrections to single-particle energies. The 
major part of particles has energies e ^ \/j,\ and is charac- 
terized by a multiple occupation N(e) ~ T/e > 1. This 
manifests itself in the dependence of interaction matrix 
elements on N(e). 

In this regime there are two energy scales (T and fi <C 
T) characterizing the distribution of particles. What are 
the particle energies that determine the many-body de- 
localization? Let us apply Eq. ^ to particles with ener- 
gies e ~ E in the energy interval of width ~ E. A typical 
value of the energy spacing is A t = l/v(E)Q(E), and 
the typical matrix element of the two-body interaction is 
enhanced due to a multiple occupation of single-particle 
states: U t = [g / ((E)](T / E) . The number Ni of occupied 
levels at distances smaller than ((E) from a given state 
is Ni(E) - Ev(E)C{E) and Eq. ([§ takes the form 

N^t/At- gTy(E)((E). 



We now use Eq. (6|) and analyze the phase diagram of 
weakly interacting disordered ID bosons. It is convenient 



Then, using the high-energy density of states v(E) — 
1/ \/Ei,( 2 E and localization length ((E) = C*E/E^ we 
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find a remarkable result: the criterion of derealization 
does not involve the single-particle energy scale! The 
transition temperature follows from the relation gT c ~ 

E*C* = \Jh 2 El/m. (The fact that this is valid for all 
energy scales suggests that the expression for the critical 
temperature/disorder can contain a prefactor logarithmic 
in 7, which we neglect). In terms of the parameters t and 
k the relation for the critical disorder becomes 

« c = t 2 /3 7 l/3. 7 -l/2 <f<7 -l (g) 

Equation (|8|) shows that k c > 1 in the entire temperature 
interval ^~/ 2 < t < 7 -1 . Hence, close to the transition 
we have E„ > ng and characteristic particle energies are 
~ > E„ > ng. This justifies our assumption that the 
major part of particles occupies high-energy states, and 
the interparticle interaction affects neither the occupa- 
tion numbers N(e) nor the chemical potential u. 

Consider now T = 0. For k ^S> 1 the boson density 
is fragmented into "lakes". Lake number i is formed by 
Ni bosons in the single-particle eigenstate \i) which is 
characterized by energy and localization length Q ss 
C*. The energy cost E^ of bringing an extra particle to 
this lake is enhanced by the repulsive interaction between 
the bosons, J3j w tj + gN/C*, and it should be equal 
to the global chemical potential /i measured from the 
lowest low-energy state e = 0. The bosons thus occupy 
only states below the chemical potential, £j < /1, with 
the occupation numbers Ni ~ (u — Ci)C,*/ g, and as long 
as u < E. t only low-energy states are occupied. The 
density of these states is = 1/E*£*, and thus the mean 
density of bosons is related to the chemical potential by 
n = fi 2 /2gE*. The chemical potential can be expressed 
in terms of the parameter k as u = E^/^/R, i.e. the 
chemical potential is indeed smaller than E* provided 
that k > 1. In this regime only a small fraction (~ u/E*) 
of low-energy states is occupied and neighboring lakes are 
separated by a distance ~ C*V^j while their size is 
~ C* ■ Using equations ([l])-(|5]) one can show that a typical 
lake (e, ~ /i) contains Ni ~ nC,*\fn — l/-v/7 bosons. 

The fact that 1(k) ^S> C* implies that for k ^ 1 the 
system is a strong insulator: the coupling between dif- 
ferent lakes is exponentially small in k. As soon as k is 
reduced to the value of the order of unity, the distance 
between neighboring lakes becomes of the order of their 
size C* and the intcrlakc coupling drives the system to the 
fluid state. So, the insulator-fluid Berezinskii-Kosterlitz- 
Thouless transition at T — occurs for k c ~ 1. 

Note that at the lower bound of the temperature in- 
terval in Eq. t ~ 7 _1 / 2 , the critical disorder is also 
K c ~ 1. Therefore, one expects that in the entire tem- 
perature range t < 7 -1 / 2 it remains k c ~ 1. 

Why the insulating state of bosons is stable at these 
temperatures as long as k > 1, i.e. E* > ng? Let the 
disorder be as str ong and reduce the temperature be- 
low y/E*Td (t < ■ Under these conditions we have 
\fjt\ < E* and only a few bosons are hosted by high-energy 
states. For high-energy bosons (e > \a\) the condition ^ 



is not satisfied because their density is too small for the 
many-body derealization. It turns out that the main 
body of the bosons, which occupy low-energy states, also 
forms an insulator. Indeed, for low-energy states we have 
e < E„ and the number of channels N\ ~ ez^(e)£(e) 
is smaller than unity, since z^(e)£(e) ~ E^ 1 . In other 
words, most of the particles occupy single-particle states 
which are separated from each other by distances exceed- 
ing the localization length £». This causes exponential 
reduction of Ut and, according to Eq. (pj, prevents der- 
ealization. Finite temperature in this situation does not 
lead to any increase in the phase volume for available 
transitions. By contrast temperature fluctuations of the 
number of particles in each lake lead to a growth of the 
energy mismatch by an amount ~ y/Tg/£*, which fur- 
ther suppresses the probablity to hybridize several states. 
The physical situation is somewhat similar to the one de- 
scribed in RefPP, where it was demonstrated numerically 
that in a finite-width band with less than one state per 
localization length the insulator remains stable with re- 
spect to the interactions at arbitrarily high temperatures. 

If the disorder is weak, k < 1, the chemical potential 
is determined by the interaction and always exceeds E*. 
As a result, we have N\{T) > 1. The condition ^ then 
indicates that the insulator is unstable, and one deals 
with the fluid state. 

The arguments presented above indicate that the phase 
transition line has to be almost horisontal in the region 
t < 1/^/7- We can also realize that at t — > the line 
K c (t) should terminate at the quantum phase transition 
point k — 1. Indeed, assuming K c (t = 0) < 1 we arrive 
at a contraidiction as the underlying superfluid phase 
at T — has delocalized excitations (phonons) at low 
energies. On the other hand, the assumption of k c > 1 
is also not consistent with the criterion (|6| as all of the 
excitations are localized and the temperature should be 
finite to provide a finite density of the excitations. 



k = E*/(gn) 



1 




t = T/(gn) 

FIG. 3: Phase diagram for weakly interacting disordered 
bosons. The green line shows the zero temperature alge- 
braic superfluid. The red part of the curve n c (t) follows from 
Eqs. Q and (JsJ> , and the blue part is an estimate for the 
low-temperature regime. 
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This completes our description of the finite temper- 
ature fluid-insulator phase transition for ID interacting 
bosons. The phase diagram is presented in Fig. [3] De- 
tailed calculations and a more accurate description of the 
temperature range t < 7 -1 / 2 (blue part of the curve K c (t) 
in Fig. [3]) will be presented elsewhere. 

IV. DYNAMICS OF EXPANSION 

The observation of the fluid-insulator phase transi- 
tion described above is feasible in experiments with cold 
bosonic atoms in the ID geometry. For electrons in 
solids one measures e.g. the DC conductivity. In quan- 
tum (neutral-atom) gases the analysis of transport prop- 
erties is based of the dynamics induced by significant 
external perturbations which may drive the system far 
from the initial state. We believe that the localization- 
delocalization transition can be identified in the dynam- 
ics of expansion of disordered bosonic clouds released 
from the superimposed trapping potential. 

In an array of harmonically trapped quasi- ID tubes 
with about a hundred of atoms per tube, the density n 
can be made ~ 10 4 cm -1 so that the length L of each tube 
is ~ 100 and the temperature of quantum degeneracy 

is of the order of tens of nanokelvins. Tuning the inter- 
action strength by Feshbach resonances or by variations 
of the tight transverse confinement one can achieve the 
interaction energy ng of the order of nanokelvins or even 
smaller and make 7 ~ 10 _1 — 10~ 3 . For the correlation 
length of the disorder a ~ 0.3 jum as in the experiment 7 
and typical values of Uo 7S , the localization length is 
in the range of 1 — 10 microns and the energy E* > ng. 
The dimensionless strength of the disorder, k §5§, ranges 
from about unity to large values and one can to study 
all temperature regimes of the phase diagram in Fig. [3] 
Note, however, that the conditions described here are 
quite different from those in the experiment!^, where the 
interaction energy was greatly exceeding E* and the sys- 
tem was deeply in the fluid state. 

Switching off abruptly the ID trap but still keeping 
the disorder (and the transverse tight confinement) like 
in the experiments^^, is expected to cause the expansion 
of the cloud. Not very far from the fluid-insulator transi- 
tion, the localization length £ of single-particle states in 
the initial cloud is much smaller than its size L. Under 
this condition, the size rapidly increases by an amount 
of ~ C "C L. If the entire initial cloud is in the insulator 
phase, i.e. locally k > tt c (t), the expansion then stops 
(see Fig. [4^,). If the central part of the cloud is in the 
fluid phase (while the outskirts are insulators) , two-body 
scattering processes induce further expansion of the cen- 
tral part. In this stage the expansion is a slow diffusive 
process governed by the diffusion equation: 

dn 9 s dn 
where r denotes the time and t is the dimension- 



less temeperature Q. The diffusion coefficient D(n,t) 
strongly depends on the local density n(x) and vanishes 
at the border of the insulating phase, i.e. for the den- 
sity n c (t) = E*/(gK c (t)). Hence, the local decrease of 
the density stops when it becomes equal to n c (t). Thus, 
the density profile of the expanded cloud represents a 
plateau with n = n c (t) in the central part and the initial 
insulating wings with n(x) < np^ (see Fig. ^3p). 

Moderately far from the transition, i.e. at (n — 
n c )/n c r-j 1, the diffusion coefficient can be estimated as 
follows. The occupation A/j of the state \i) changes by ~ 1 
on a time scale of the order of the inverse matrix clement 
fi/Ut ~ fi^/gNi. The distance of a typical hop is ~ £. 
Therefore, on a time scale tq ~ Nih/Ut ~ fr£/<? each bo- 
son will move by a distance ~ Q, and D ~ C 2 / T o ~ Cfi 1 /^- 
For high-energy states the localization length at typi- 
cal energies e is £(e) = ^e/E*, and the diffusion coef- 
ficient is given by D ~ e(*g/TiE*. Locally, the decrease 
of the density ceases when it reaches the critical value 
n c (t). Thus, for the initial central density comparable 
with n c {t), the evolution of the cloud to the final shape 
of the plateau and wings (see Fig. B| requires a charac- 
teristic time ~ L 2 /D. Since L ~ (e/mw 2 ) 1 / 2 , where w 
is the initial trap frequency, we obtain 

1 hE, 1 (ng\ k 3 / 2 
T * 7 TH- ( 10 

However, in the temperature interval < T < T4 

most of the particles, which have energies e ~ \fj,\ = 
T 2 /Td, and particles with e ~ T will expand with dif- 
ferent velocities. So, the redistribution of particles in the 
course of the expansion may become important and it 
can slightly modify the above estimate for the time r*. 

The time by far exceeds the time u^ 1 of ballistic 
expansion of the clean thermal cloud. For example, near 
the lower bound of Eq. Q, i.e. for T ~ T d and K ~ 7~ 1//3 , 
we have ~ Td/hto 2 . For typical values Td ~ 50 nK and 
u) ~ 10 Hz this estimate yields t, ~ 1 s. Moving to the 
lower bound of Eq. d8| where T ~ Td^/j and k ~ 1, the 
time r* reduces by a factor of »/7, i.e. it is of the order 
of 0.1 s. As we see, these time scales are such that one 
can observe the evolution of the expanding cloud from 
the very beginning until the final density distribution is 
formed. The expansion of the fluid part of the cloud will 
also be slow in the low-temperature regime, T < Tdy/^y 
(t < 7 -1 / 2 ). The analysis of the diffusion coefficient for 
this case is beyond the scope of the present paper. 



V. PHASE DAGRAM IN HIGHER 
DIMENSIONS 

We conclude our discussion of the properties of disor- 
dered interacting bosons with a brief sketch of the phase 
diagram in higher dimensions. As is well known, in the 
absence of disorder bosons form a superfluid (algebraic 
superfluid in 2D) below a critical temperature T c . At 
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FIG. 4: Manifestation of the many-body localization- 
delocalization transition in the expansion of a quasi-lD cloud. 
In (a) the the entire initial cloud is an insulator, whereas in 
(b) initially the central part is in the fluid and the outskirts 
in the insulating phase. Initial and final shapes of the cloud 
are shown in black and red, respectively. 
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FIG. 5: Phase diagram for two-dimensional weakly inter- 
acting bosons. The black curve shows the thermodynamic 
Berezinskii-Kosterlits-Thouless transition. The red curve is 
the many-body localization-delocalization transition. The 
blue dashed curves indicate phase boundaries which can be 
ruled out on rather general grounds (see text). 



high temperatures, T > T c , the clean system is a nor- 
mal fluid. The superfluid survives a sufficiently weak 
disorder, but the superfluid transition temperature T c de- 
creases with increasing the strength of the disorder and 
vanishes at a critical strengtiP&l (see the black curve and 
brown point in Fig. [5j. What is the state of the dis- 
ordered system at T > T c , i.e. above the black curve 
in Fig. [5] ? It follows from our previous discussion that 
bosons can form either the normal fluid or the insulat- 
ing state. The suggested phase diagram is presented in 
Fig. [5] 

This sketch can be justified in the following way. First, 
it is safe to assume that at sufficiently strong disorder 
bosons are in the insulating state. The second obser- 
vation is that at the critical disorder (brown point in 
Fig. |5| and arbitrarily small but finite T one should ex- 
pect the normal fluid rather than the insulator. Indeed, 
the zero-temperature insulator can be thought of as a 
system of superfluid lakes which are separated from each 
other and have uncorrelated phases. A typical size of 
such a lake increases with decreasing the strength of the 
disorder and diverges at the critical strength. It means 
that although excitations in the insulating phase at T = 
are always localized, their localization length can be ar- 
bitrarily large. Therefore, at any finite temperature and 
interaction strength there is a vicinity of the critical dis- 
order, where the insulator is unstable with respect to the 
many-body derealization. As a result, the insulator - 
normal fluid phase boundary can not follow the lower 
blue curve on Fig. [5] 

On the other hand, the normal fluid can not be sta- 
ble at T = . Indeed, the wavefunctions of low-energy 
single-particle states have to be localized, otherwise par- 
ticles can not avoid Bose-Einstein condensation and the 
system would become superfluid. At extremely low but 
finite temperatures, the density of thermal excitations is 
vanishingly small and the interaction between them is un- 
able to delocalize many-body wavefunctions. This rules 
out the phase boundary following the upper blue curve 
in Fig. [5j Thus, the only possible option is represented 



by the red solid curved 

The arguments given above apply to both 2D and 3D 
cases. However there is a big difference, since in 2D all 
single-particle states arc known to be localized, and the 
strong-disorder phase is a true insulator. This means that 
the diffusion constant is zero even at finite temperatures. 
At the same time, in three dimensions high-energy states 
can be extended. As soon as the extended states appear 
they can host thermal excitations and thus allow expo- 
nentially small but finite diffusion. From this point of 
view, there is no qualitative difference between the insu- 
lator and normal fluid, and the red curve in 3D represents 
a crossover rather than a true phase transition. 



VI. CONCLUDING REMARKS 

A remarkable possibility to compare disordered inter- 
acting bosons with non-interacting ones is offered by a 
7 Li atomic gas, where the coupling strength g can be 
varied by a Feshbach resonance from practically zero to 
large ppositive values^. In the ID case, achieving the 
strongly interacting regime where 7 > 1, is expected to 
present new transparent physics. For 7 — > 00 as in the 
recent cesium experiments, the bosons become impen- 
etrable and show a clear analogy with non-interacting 
fermions. All single-particle states are then localized irre- 
spective of the strength of the disorder. At intermediate 
values of 7 one expects a peculiar interplay between the 
interparticle interaction and temperature. For T 3> ng 
the situation should be the same as desribed in Section 
III for the high-temperature regime, T T<j, and in this 
sense equation ^ is universal. In the other extreme, 
T < ng, an increase in the interaction strength first leads 
to the localization-delocalization transition at small 7, 
but then causes a reentrance to the insulating phase at a 
critical disorder-dependent value of 7. This behavior is 
expected because at T — even an infinitesimally small 
disorder leads to the appearance of an insulating (Bose 
glass) phase if 7 is sufficiently large^S. The comparison of 
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the insulating phases emerging at small and large 7 with 
each other is supposed to shed new light on the structure 
of Bose glasses. 
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